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We discuss four general features of force-free evolution: (1) The spatial spread of any packet 
changes with time in a very simple way. (2) Over sufficiently short periods of time (whose duration 
is related to the spread in momentum of the packet) the probability distribution moves but there 
is little change in shape. (3) After a sufficiently long period (related to the initial spatial spread) 
the packet settles into a simple form simply related to the momentum distribution in the packet. 
In this asymptotic regime, the shape of the probability distribution no longer changes except for 
its scale, which increases linearly with the time. (4) There is an infinite denumerable set of simple 
wave packets (the Hermite-Gauss packets) that do not change shape as they evolve. 



I. INTRODUCTION 

The behavior of free wave packets as they evolve in 
time is striking, even alarming, to many students of quan- 
tum mechanics. A wave packet usually changes shape 
and always eventually spreads out without limit. How 
different this is from the behavior of a free classical parti- 
cle! Reconciliation of the classical and quantum descrip- 
tions is part of the bigger question of the interpretation 
of quantum mechanics, but the dynamics of free wave 
packets is still of considerable interest in introductory 
quantum mechanics. 

Most texts derive the propagator 

K(x,x',t) = yjm/2mht exp[im(a; - x') 2 /2M], (1) 

such that the evolution of an arbitrary initial wave packet 
ip{x,0) is exactly 

/oo 
K(x ) x',t)ip(x',0)dx'. (2) 
-oo 

Unfortunately, there are not many wave packets for which 
this integral can be easily evaluated and the only case 
commonly treated is the Gaussian packet. Therefore the 
evolution of free wave packets remains obscure to many 
students. 

Ehrenfest's result that (p) is constant and that (x) — 
(x) +(p)t/m is well known. There are four other insights 
into the free evolution of wave packets that are simply 
derived and greatly improve the general understanding 
of this evolution: 

1. The spread of any wave packet, as measured by 
A x = ((x — (i)) 2 ) 1 / 2 , varies in time as 

A x = v /A 2 im + (<-t mm ) 2 A 2 /m 2 , (3) 

where A m i n is the minimum spread and A p = ({p — 
(p)) 2 ) 1 / 2 does not change with time. If the wave function 
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is real everywhere (which can be only for an instant) then 
the packet has its minimum spread A m i n at that instant. 

2. Over sufficiently short periods of time, much less 
than mh/A 2 , the wave packet moves (with speed (p)/m) 
without much change in the shape of the probability dis- 
tribution. 

3. The asymptotic evolution, for times \t\ 3> mA 2 /h, 
is to a form simply related to the momentum distribution 
in the packet. This means that for any initial packet for 
which the momentum distribution (the Fourier transform 
of the wave function) can be easily evaluated, the evolu- 
tion of the packet in this asymptotic period is accessible. 

Note that these two time regimes can not overlap, be- 
cause of Hciscnbcrg's uncertainty relation A^Ap ^ h/2. 

4. There is an infinite number of simple wave pack- 
ets that do not change shape (apart from the inevitable 
spreading) as they evolve. These wave functions have the 
form of a Gaussian multiplied by a polynomial in x. 

In a recent note [I], it is claimed that "a free particle 
wave packet of any shape becomes approximately Gaus- 
sian after a period of time" . This claim is false. An obvi- 
ous counter-example is provided by any odd wave packet: 
it will remain odd as it evolves freely and therefore can 
never approximate a Gaussian. Furthermore, the claim 
is not consistent with Properties 3 and 4 above. This 
note is discussed further in the appendix. Here we will 
establish the four properties and give some examples. 

We use the term wave packet to indicate that the wave 
function is sufficiently confined in space for (x 2 ) to exist. 

II. THE EVOLUTION OVER SHORT PERIODS 

Starting with the well-known formula[2] for the evolu- 
tion in terms of an initial momentum distribution <fi(p), 

1 f°° i 1 
i/)(x,t) = -== exp[-(px- — p 2 t)](j)(p) dp, (4) 
\Z2irh J~oc n, 2m 

expand 3 p 2 about the mean value p = (p). Thus 
p 2 = p 2 + 2pp + (p — p) 2 . We expect <f>{p) to be small 
if (p - p) 2 > A 2 and therefore, if A 2 p t/2mh < 1, the 
term in (p — p) 2 will make a negligible change to the 
phase of the exponential. In other words, the integrand 



2 



contains a factor exp[t(p — p) 2 t/2mh] that will be very 
close to unity throughout the region where (p(p) has any 
significant magnitude. Hence 



ip(x,t) 




! XP[-P(X 



Pt 



(5) 



The wave function translates, with speed p/m, and no 
change in shape, just a change in phase. A more rigorous 
approach to the error in this approximation is given in 
Section VII. 

Here we have taken the initial time to be t = 0, but 
this analysis can be applied at any time. Since A p does 
not change, the time scale mh/2A 2 for changes in shape 
applies throughout the evolution. 

If the wave function has any discontinuity then A p = 
oo and there will be very rapid changes. An example of 
this behavior will be considered in Section VI. 



III. THE EVOLUTION FOR LARGE TIMES 

For the propagator in Eq. 0, write the term (x — 
x') 2 in the exponent as x 2 — x 2 — 2(x — x)x' + (x' — x) 2 , 
with x = (x)t=o- We expect ip(x' , 0) to be small if (x' — 
x) 2 3> A 2 and therefore, for all times (future or past) 
with \t\ S> mA 2 /h, the term in (x' — x) 2 will make a 
negligible change to the phase of the exponential. Hence 



ip(x,t) 



^ — tt exp — —lx — x ) 



exp [ — j~^( x ~ x)x'] ip(x' , 0) dx' 



j t exp[^(x 2 -x 2 )]0(m(x-x)A), (6) 



where 



Hp) 



V2irh J- 



exp(—ipx'/h)ip(x',0)dx', (7) 



which is the initial momentum wave function of the 
packet. A more rigorous derivation of the error in 
this approximation is given in Section VII. Historically, 
this asymptotic form was important in the theory of 
scattering. [5] 

The asymptotic probability distribution is therefore 
(m/t) \<p(m(x — x)/t)\ 2 . This has a simple interpretation: 
m(x — x)/t is the momentum required for the particle to 
be at position x at time t, if it left x at time t — 0. So the 
probability of finding the particle at x is proportional to 
the probability that it had the right momentum to get 
there (ignoring the fact that the initial distribution was 
spread over a distance of order A x about x). 

Introduce the times t p = mh/2A 2 and t x — 2mA 2 jh. 
Then t p ^ t x (from the uncertainty relation) and equal- 
ity implies a Gaussian packet. In general, t p gives the 



time scale for changes in shape (including the changes 
in spatial scale), and t x gives the time scale to move to 
the asymptotic period where there is no further change 
in shape but \ip\ expands uniformly with time. 



IV. THE EVOLUTION OF THE SPREAD A x 

The fact that every free packet spreads exactly as in 
Eq. (J3) has appeared in at least one textbook [5], but 
the derivation is simpler if the concept of the total time 
derivative of an operator is used[B]. The total time 
derivative, d t A of any operator A is defined to be 



d t A = d t A + (i/h)[H,A], 



(8) 



where H is the Hamiltonian operator. For the free 
particle this leads to the operator equations of motion 
mdtX = p and dtp — 0. From the general property 
that d t (A) = (d t A), it follows that md t X = P and 
d t P = 0, where X = x — (x) and P = p — (p). Fur- 
thermore, md t X 2 = 2R, where R = (PX + XP)/2, and 
mdtR = 1 ■ Since (P 2 ) is constant, it follows that 



m(R) = (P 2 )t + m{R) t=0 = (P 2 )(t - t mm ), (9) 



where t r 



md t (X 2 ) = 2(R) to give 



-m(R)t=o/(P 2 )- Now we can integrate 



(x 2 



l (P 2 )(t-t mm ) 2 + (X 2 }, 



(10) 



which is equivalent to Eq. ([3]). Thus, for any wave packet, 
one can calculate the time of minimum spread from (P 2 ) 
and (R)t=o- A wave packet has its minimum spread at 
the time when (R) = 0. 

If the wave function tp is real at some time, the i in 
Schrodinger's equation dtip = —t(h/2m)d 2 ip shows that 
ip will immediately become complex. Also (p) = be- 
cause (p) = fippipdx — f{pip)*ip dx = —j(pip)ipdx — 

— (p). Similarly, (px) = Jippxipdx — J (pip)* xtp dx = 

— / (pip) (xtp) dx——(xp). That is, (R) = and therefore 
the wave packet has its minimum spread at that instant. 
This also shows that a free wave packet can be real only 
for one instant. In fact, the probability is symmetric 
in time about that instant, because if ip(x, i) satisfies 
Schrodinger's equation then so does ip*(x,2to — t) and 
since these two functions are equal at t = to if ip(x, t ) is 
real, they must be equal at all times. 

Given a wave packet £ that is at rest (in the sense 
that (p)^ = and therefore (x)j is constant) we can 
make a moving packet by multiplying the wave func- 
tion by exp(ipx/h). With ip(x,to) = exp(tpx/h)!;(x,to), 
it follows that pip = pip + exp(zpx/h)p^ and therefore 
(p)tl) — P- Similarly, (p 2 )^ = p 2 + (p 2 )^ and hence 
(Ap),/, = (Ap)j. Also, p(xip) — pxip + exp(ipx/H)p(xt;) 
and therefore (px + xp}^ — 2p(x}^ + (pi + xp)^. Thus, 
(R)^ = (R)^ and it follows that the evolution of the 
spatial spread is the same as for £. The evolution of ip 
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can be obtained from that of £ by applying a Galilean 
transformation [7] : 

TP(x,t)=exp(^(px-p 2t -^))£(x- p t -^,t). (11) 

It is easy to check that this ijj satisfies Schrodinger's equa- 
tion if £ does, and that ip(x, to) = exp(ipx/h) to). 



V. THE HERMITE-GAUSS WAVE PACKETS 

A vast array of time-dependent solutions of 
Schrodinger's equation (SE) can be found through 
the use of invariant operators [5] . For any operator A, 

(H - ihd t )Atp = ([H,A] - ihd t A)ip = -ih(d t A)4>. (12) 

If dtA = at all times, we say that A is invariant. If 
an invariant operator is applied to any solution of SE, 
then from Eq. ( fl2| the result will also satisfy SE. Af- 
ter being multiplied by a suitable time-dependent factor 
(if necessary), cigenfunctions of invariant operators will 
satisfy SE and the eigenvalue will be constant. [If the 
eigenfunctions are degenerate, appropriate linear com- 
binations, with time-dependent coefficients, may be re- 
quired.] 

In the present case of a free particle, a useful invariant 
operator is 



b = 



p{t — it), t real 



(13) 



The eigenfunction of b with eigenvalue zero (so that (x) — 
(p) = 0) is easily found using p = —ihd x : 



X(M) = (-) 



1/4 _ 



: exp(- 



V(t - it) * x 2h{t-iTy 



(14) 



This is Schrodinger's well-known free Gaussian wave 
packet. The factor (t — it)^ 1 / 2 does not come from the 
eigenvalue equation b\ = 0, but it is clear from the prop- 
agator in Eq. ([I]) that this factor is required to satisfy SE. 
Propagators always satisfy SE and with x' — the prop- 
agator in Eq. ([7J) has the same form as x with r = 0. 
Furthermore, any solution of SE is still a solution if a 
constant (such as —it) is added to t; thus \ satisfies SE. 

Applying b' = mx — p(t + it) to the eigenfunction x 
will give another solution to SE (because is also in- 
variant) , and clearly this must have the form of the same 
Gaussian multiplied by x. Applying ¥ again gives a so- 
lution that has the form of the same Gaussian multiplied 
by a quadratic polynomial in x, and so on. We will now 
show that the sequence of polynomials generated in this 
way are Hermite polynomials. 

First write exp[tmx 2 / 2h(t — it)] in Eq. ( 14 1 as 
cxp[imx 2 /2h(t + it)] cxp[—m,TX 2 /h(t 2 + t 2 )] and then, 
since b' = mx + ih(t + it)8 x , 



b' exp[ 



1 2 

^imx 
h(t + it) 



] f = ih(t + it) exp[ 



h(t + it) 



d x f (15) 



for any function f[x,t). Now we can apply this result 
repeatedly, starting with x, to give 



Xn(x,t) (X 



>vmx 



{t + lT) n 

y/(t-lT) ^ l h{t + lT) 



exp[ 



d: exp(--), (16) 
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where we have introduced the time-varying length scale 
7 = [h(t 2 + t 2 ) /mr] 1 / 2 . The generating relation (Ro- 
drigue's formula) for the Hermite polynomials is 

(d/dZ) n exp(-a = (-l) n H n (0 exp(-a, (17) 
and using this in Eq. (fl6| gives 



(t+iT) n j n ^imx 2 _ iB 2/~2 x 
Xn{x,t) oc y - 7 7 -— — exp 2 e n . 18 

y/(t - IT) h{t+lT) 7 

Separating the real and imaginary parts of the first ex- 
ponent, the normalized sequence of solutions to SE is 



Xn(x,t) = 



7T-V 4 

V7 2 n n 



e *(fl+(n-l)/J) e - 2 /27 2 Fn (^) ) ( 19 ) 



where 9 = tx 2 /2Tj 2 and e 4/3 = [{t - vr)/(t + it)] 1 / 2 . The 
probability distributions of the first few are: 



lx(M)| 2 = 
lxiOM)l 2 = 
lx 2 (M)| 2 - 



l 



7V7T 

2 x 2 



cxp(— x 2 /j 2 ) 



j^/ir 7 



exp(— x /7 ) 



(20) 
(21) 



2j^/tt 



(l-2x 2 /7 2 ) 2 exp(-x 2 /7 2 ).(22) 



These do not change their shape as they evolve; only 
the scale of x changes with time. It would serve little 
purpose to show a graph of these probability distributions 
because they have the same spatial form as do the energy 
eigenfunctions of the harmonic oscillator. The relation to 
the harmonic oscillator is made clearer in Section VII. 

The fact that the Hermite-Gauss packets do not change 
shape is consistent with the asymptotic results in Section 
III because the Fourier transform of a Hcrmite-Gaussian 
is also a Hermite-Gaussian. 

1 Z" 00 

= e~ z / 2 H n (z)e- ipz dz = i- n e- p ' 2 H n {p). (23) 

2TT J-oo 



VI. EXAMPLES 

Our first example is a wave packet which is smooth 
(infinitely diffcrcntiablc) everywhere and goes to zero ex- 
ponentially at large distances, but does change shape as 
it evolves. We use the second derivative of Schrodinger's 
packet x m Eq.(14l. Thus X2(x,t) cx 9 2 x(a;,^), which 



also must satisfy SE. It has the form 



X2{x,t) = Nk 6 (2k 2 x 2 - l)e" K 



(24) 
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FIG. 1: The probability distribution \ip(x, t)\ 2 for the wave 
packet in Eq.(24l, whose initial probability distribution is 




shown by the solid curve. Only half of the packet is shown, 
because it remains symmetric. The probability is shown 
in the period before the asymptotic regime for the times 
t = 0,0.1,0.2,0.4,0.6,0.8,1.0 in units of r. The distance x 
is in units of (ftr/m) 1//2 . 
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FIG. 2: The probability distribution for the same wave packet 
as m Fig. [l] for some times (t = 3, 4, 6, 16) when the packet is 
approaching its asymptotic form. In the asy mpt otic region, 
t\tp(x,t)\ 2 is a function of x/t only, as in Eq.(25l, so we use 
x/t as the horizontal variable in this figure, which hides the 
substantial spreading of these packets. For all t > 7t x w 16r 
the exact distribution is indistinguishable from the asymp- 
totic form; so the graph for t = 16 stands for all later times. 



where k 2 — —im/2H(t — it) and N is the unimportant 
normalizing constant (3v / 7r/32) _ 2 . Initially, the wave 
function is real. It is also symmetric about x = and 
will maintain that symmetry. As shown in Fig. [I] there 
is a hump centered on the origin with a lower hump on 
either side. We will see that as it evolves, the central 
hump diminishes and eventually disappears while the 
outer humps move further out and spread. The initial 
packet has A 2 = Ihr/Qm and A 2 = 5/mti/2t, and hence 
t p — t/5 and t x — 7t/3. We expect that the shape (of 
the probability distribution) will not change significantly 
over periods much less than t p . Initially, no change can 
be seen in the graph until t ss t p /3 ~ 0.06r. The maxi- 
mum rate of change appears to be at about t = 0.6 and 
x = when visible change starts for St w t p /30. 

Taking the Fourier transform of \2 leads to the asymp- 
totic form 



X2OM) ~ N '^ ex P ' 



rnrx 
2ht 2 



(25) 



with TV' = 2(m 5 T 5 /97r/i 5 ) 1 / 4 , which we expect to be valid 
for t > 4. No error in this approximation can be seen 
in Fig. [2] after tm7t x . 

As a second example, consider the square wave packet, 
with 



I 0, \x\ > a/2. 



(26) 



The momentum wave function is, from Eq. ([7]), 

m-- 



a sm(ap/2K) 



2irh ap/2h 



(27) 



and therefore, for t 3> ma 2 /h, 



2-Kiht 



exp[ 



.irax 2 sm.(amx /2M) 



2ht 



amx/2M 



(28) 



The exact evolution can be carried out, using Eq. (§, 
in terms of error functions with complex argument, or 
Fresnel functions with real argument. Fig. [3] shows the 
exact probability distribution for some times well before 
it reaches the asymptotic realm. Fig. [4] shows the exact 
probability distribution for some times when the exact 
distribution is approaching its asymptotic form. 

In this case the initial value of A 2 is a 2 /12 so that t x = 
ma 2 /6h. Fig.|4]shows that the asymptotic form is a close 
approximation for t > 3t x . The short-period theory in 
Section II does not apply because A p = 00, and there are 
rapid changes during the early evolution. Discontinuous 
wave functions are unphysical; an infinite energy would 
be required to produce them. This does not mean that 
they cannot be useful as mathematically simple states 
that can be closely approximated by physical ones. 



5 



0.5 



.0.1 



\0.01 4 

i \ 0.001 



.! ' If 



2 \ 



0.0 0.2 0.4 0.6 0.8 



1.0 



FIG. 3: The probability distribution \ tp(x, i)| 2 for the initially 
square wave packet in Eq.(26l, for several values of the time 
(t — 0,0.001,0.01,0.1) in the early period when the packet 
shows complicated behavior, but does not spread appreciably. 
The distance x is in units of a, the width of the initial packet. 
The times are in units of ma 2 /H and are indicated by the 
numbers against the graphs. 



Now make use of the Cauchy-Schwarz inequality 
I Jf(p)g(p)dp\ 2 < J\f(p)\ 2 dpJ\g(p)\ 2 dp, with g(p) = 
(p-p)Hp) an d \f(p)\ = I sm[(p - p) 2 t/4mh]/(p - p) |, 
and the integral sin 2 (z 2 ) / ' z 2 dz = y/jr, to show that 



\Sip\ 2 < y/t/irmhZA 2 



(30) 



This is a rigorous bound on the change in the wave func- 
tion (apart from the shift pt/m). For \Sip\ to be small 
compared with ip with I?/)) 2 ~ l/2A x (from the normal- 
ization of ip), we require tA^/irmfv' <C 1/4A 2 and hence 
tA 2 /Trmh 3 <C 1/4A 2 A 2 < l/h 2 (from the uncertainty re- 
lation). That is, we require t <C irmh/A 2 , in agreement 
with Section II. 

(b) For a more rigorous approach to the error in 
the asymptotic approximation, write the exact form of 
Eq. ([6) as 



ip(x,t) 



where 



vf exp[ 



^-y-x 2 ))^{x-x))+5^ (31) 
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FIG. 4: The probability distribution for the initially square 
packet for some times (t = 0.1,0.2,0.5) when the packet is 
approaching its asymptotic form. We show t\ip{x, t)\ 2 as a 
function of x/t. For all t ^ 0.5 the exact distribution is in- 
distinguishable from the asymptotic form; so the graph for 
t — 0.5 stands for all later times. 



VII. FURTHER CONSIDERATIONS 

The material above is simple enough and sufficiently 
complete to be presented in an elementary course in 
quantum mechanics, but some related topics may be 
helpful and may serve as extension material for students. 

(a) For a more rigorous approach to the short-time 
behavior, write the exact form of Eq. ^ as ip(x, t) = 
exp(—tp 2 t/2mh) ip(x — pt/m, 0) + dip, where 



Sip(x, t) 



exp( W x 



exp[-p(x 



— )] (cxp[ ~ l( - P ~ P ^ ] 
m 2mh 



i)(f>(p)dp. 



Sip = 



2-Kiht 



im . 2 

exp W* 



X 2 )] 



f(x')g(x')dx', (32) 



with 



f(x') = exp[^(x 



x)x 

g(x') = (x'-x)iP(x',Q). 



Now the Cauchy-Schwarz inequality gives 
\Sip\ 2 ^ yfmF/^W^A 2 x . 



(33) 
(34) 

(35) 



This is a rigorous bound to the error in the asymp- 
totic approximation. For \8ip\ to be small compared 
with ip we require |5V>| 2 "C 1/2A X , which leads to t ^> 
(4/7r) 1 / 3 mA 2 , in agreement with Section II. 

If the wave packet has discontinuities, as in the square 
wave in Eq. ( 26 1 , this derivation of the error is applicable 



only when the initial time is taken to be the instant when 
the discontinuities exist. At any other time, A x does not 
exist (because m 2 d 2 (X 2 ) = 2(P 2 ). as in Section IV, and 
(P 2 ) = oo.) The divergence of (x 2 ) can easily be seen 
for the asymptotic wave function of the initially square 



packet from Eq. ( 28 1 



(c) The relation of the Hermite-Gauss solutions for free 
particle motion to the energy eigenfunctions of the har- 
monic oscillator can be further elucidated as follows. Sep- 
arating the real and imaginary parts of the exponent in 
X in Eq. ([T4| gives 



(mr /nh) 1 ^ 4 , itx 
X{x,t) = y ' ' exp( 



2ry 



x 2 



(36) 



(29) w ith 7 = [h(t 2 + t 2 ) /mr] 1 / 2 . Similarly to the derivation 



of Eq. (15 1, one can show that 



•be 1 



i\j2hmT e %f3 a, 



(37) 
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where a = (x/j+vyp/ti) /y/2 and, as before, 9 = tx 2 /2tj 2 
and (3 = ((t - vr)/(t + it)) 1 / 2 . Then [a, fit] = 1 and 
a is the type of operator familiar from the theory of 
the harmonic oscillator. Schrodingcr's packet Xi with 
bx = 0, corresponds to the ground state of a, with 
aexp(— x 2 /2"f 2 ) — 0. Then the sequence of states ob- 
tained by successive application of a) to exp(— x 2 /2 r y 2 ) is 
the sequence of oscillator states exp(— x 2 /2j 2 )H n (x/"f) 
and corresponds to the sequence of solutions of the free 
SE obtained by successively applying b' to x- 

The Hermite-Gauss solutions can be readily general- 
ized to have non-zero values of (x) and (p) . Inste ad o f 
bx = leading to the solution for x{ x ,t) m Eq. (14 1, 
start with bx = ftXi where (3 = m(x) — (p)(t — it), and 
then repeatedly apply bK 

It is not surprising that this Hermite-Gauss sequence 
of solutions exists for the free particle; a similar se- 
quence can be found for any Hamiltonian that is at most 
quadratic in p and x, even when the coefficients of the 
terms in the Hamiltonian vary with the time. [5] 

(d) The spatial derivative of any solution of the free 
SE is also a solution. Taking successive derivatives of 
Schrodinger's packet \ in Eq. ( 14 1 gives a sequence of 



packets that also have the form of a Hermite polynomial 
multiplied by the original Gaussian and by a phase factor, 
but in this case the probability distributions do change 
shape as they evolve. 



From Rodrigue's Eq. (17), and defining Xn — d x x 



cx K n+1 e~ K x ' 



H n (nx), 



(38) 



where n — [m/2ih(t— it)] 1 / 2 . But the Fourier transform 
of this combination of a Gaussian and an exponential is 
different to that in Eq. p3j) : 



2 H n (z)e- ipz dz= t^e- p2 / 4 p n . (39) 



The momentum distribution of Xn is not a Hermite- 
Gaussian and therefore these packets change shape as 
they evolve. The first example in Section VII illustrates 
the case n = 2. 



VIII. CONCLUSION 

To see the big picture of the free evolution of a wave 
packet, it helps to extend the evolution back to the dis- 
tant past as well as forward. The centroid (x) of the 
packet moves with constant speed {p)/m and we know 
from Section IV that every packet will contract and 
then expand, so that A x will take its minimum value 
A mm at some time to- Far from the minimum, for 
\t — to\ th = mA m i n / A p , A x will vary linearly with 
time, from Eq. We also know, from Section III, that 
for \t — to I ^ tx — ^• m ^ 2 minl^ 1 the shape of the wave 
packet will not change, but the scale will change linearly 
with time (consistent with the linear change in A x ). Note 
that t h = (txtp) 1 / 2 , where t p = mh/2A 2 . Thus t h is the 
harmonic mean of t x and t p and t p si ^ t x . 

Therefore, changes in shape (other than in scale) can 
occur only in the waist region around to and such changes 
cannot be rapid over times of the order of t p , if A p exists. 
But even in this waist region, there need not be changes 
in shape, as exemplified by the Hermite-Gauss packets. 



APPENDIX: COMMENT ON MITA'S PAPER 

The only example considered in detail by MitapQ is 
the square wave packet, as in Eq. ( p6| . We have shown 
that in this case, evolves to approach the form 

sm.(/3x/t)/(f3x/t) where j3 — am/2h, and this packet is 
not Gaussian. All that has been achieved in Mita's note is 
to impose a Gaussian fit to an arbitrary wave packet [by 
making a quadratic approximation to the logarithm in 
Mita's Eq.(14)]; if the packet is nothing like a Gaussian, 
the fit will be poor. 

For Mita's parameters m = h = a = 1, the exact distri- 



bution differs from the approximate one [from Eq. (28 1] 
by a negligible amount for all t > 0.5. Mita's Fig. 3(c) 
shows only that, if you don't look too closely, (sinx/x) 2 
is not very different from a Gaussian shape; if the cc-axis 
were extended further, the next maximum in the exact 
distribution would be revealed - it does not decrease in 
relative size for later times. 
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